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Abstract. In minimal supersymmetric seesaw models, the contribution to lepton flavour 
violation from ^-penguins is usually negligible. In this study, we consider the supersymmetric 
inverse seesaw and show that, in this case, the Z-penguin contribution dominates in several 
lepton flavour violating observables due to the low scale of the inverse seesaw mechanism. 
Among the observables considered, we flnd that the most constraining one is the /i-e conversion 
rate which is already restricting the otherwise allowed parameter space of the model. Moreover, 
in this framework, the Z-penguins exhibit a non-decoupling behaviour, which has previously 
been noticed in lepton flavour violating Higgs decays. 



While the basic ingredients of neutral lepton flavour violation are close to being fully 
determined, charged lepton flavour violation (cLFV) remains to be observed despite impressive 
experimental efforts. In fact, cLVF is allowed in minimal extensions of the Standard Model but 
the associated observables are highly suppressed due to a Glashow-Iliopoulos-Maiani (GIM) 
mechanism taking place in the lepton sector. However, when these extensions are embedded 
in a larger framework like supersymmetry (SUSY), large contributions to cLFV are expected. 
The supersymmetric inverse seesaw is a very attractive framework to consider since all the new 
Physics it introduces does not exceed the TeV scale. As a consequence, it could be tested at high- 
energy through the direct production of the new heavy particles but also at the high-intensity 
frontier. Indeed, in the inverse seesaw pQ, the size of the effective operator at the origin of the 
neutrino masses is not correlated to the size of the one generating cLFV. This suggests that 
there could be a large enhancement to lepton flavour violating observables, which has recently 
been discussed in several studies [IHS]. As recently pointed out, Higgs-mediated observables can 
be signiflcantly enhanced by the presence of relatively light right-handed neutrinos [9]. Contrary 
to the Minimal Supersymmetric Standard Model (MSSM) where Z-penguin contributions are 
suppressed due to several cancellations, in the supersymmetric inverse seesaw, Z-mediated 
contributions are enhanced and even dominate the ones arising from Higgs-penguin. Indeed, 
in most of the cases, the Z-penguins provide the dominant contribution and can also exhibit a 
surprising non-decoupling behaviour. 

After defining the model, we will discuss the enhancement of the Z-mediated contribution, 
giving an approximate formula for the Z-l^-l^ vertex. We will then present our numerical results 
and conclusions. 



1. The supersymmetric inverse seesaw 

The supersymmetric inverse seesaw consists of the MSSM extended by three pairs of gauge singlet 
superfields, uf and Xi [i = 1,2,3), with opposite lepton number, —1 and +1, respectively. The 
superpotential of this model is given by 
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+ Y^'vmH^ - fiH^H^ + MR^MfX, + -fsx^.X^Xj , (1) 

where i,j = 1,2,3 are generation indices. In the above, and Hu are the down- and up- 
type Higgs superfields, Li denotes the SU(2) doublet lepton superfields. The only term that 
violates lepton number conservation is the Majorana mass term fiXijXiXj. Thus, the limit 
fJ'Xij — ^ increases the symmetry of the model, making the smallness of ^Xij natural. The 
soft SUSY breaking Lagrangian can be found in [ID]. Assuming a flavour-blind mechanism of 
supersymmetry breaking by considering universal boundary conditions at the grand unification 
(GUT) scale, like in case of the constrained MSSM (cMSSM), the supersymmetric parameters 
are 

m^ = mo, Mgaugino = Mi/2,Ai = AoI, B^^ = Bm^ = BqI. (2) 
In this framework, the 9x9 neutrino mass matrix Ai in the basis {v, z^'^, X} is given by 
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M = \ mo Mr \ , (3) 

Ml 

where mj) = -^YyVu and M/j, ^x are 3x3 matrices in family space. Assuming mz), ^x ^ Mr, 
the diagonalization leads to an effective Majorana mass matrix for the light neutrinos [llj . 
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= mlMl~\xM^^mD = ^yJ (mI)-^ fixM^^Y, . (4) 

As can be seen from the above equation, the smallness of the light neutrino masses can be 
explained by the smallness of ^x- This decorrelates the effective operator generating neutrino 
masses from the one generating cLFV since the latter does not depend on fix- Moreover having 
Hx at the eV-scale allows the Yukawa couplings to have natural values (Yi, ~ C(l)) and, at 
the same time, right-handed neutrinos to have masses at the TeV. Having light right-handed 
neutrinos and large Yukawa couplings leads, in turn, to large contributions to cLFV observables. 
By defining 

M-i = (M^)-i /XX , (5) 

the light neutrino mass matrix can be cast in a form that strongly resembles the one for the 
type-I seesaw 

= "^YjKr^Y, . (6) 

Without loss of generality, we work in a basis where Mr is diagonal at the SUSY scale. 
Moreover, in the numerical analysis, we will work in the limit where ^.x is diagonal. This is a 
conservative assumption, since, in that case, the experimental limits will constrain the size of 
the Yukawa couplings, which in turn will limit the branching ratios and decay rates of cLFV 
observables. The constraints from the neutrino data that we considered in this study can be 
found in [TOl. 



2. Z-mediated charged lepton flavour violation 

Many studies have been devoted to cLFV in the MSSM and many of its extensions. It was 
found that the 7-penguin gives the dominant contribution to 3-body cLFV decays [12^13). the 
only exception being the large tan/3 (and low mA) regime where large Higgs contributions are 
expected |14) . The same observation was recently made in the super symmetric inverse seesaw [9]. 
When the mass scaling of the different contributions is considered, the 7-penguin contribution 
to ii — 7> 3ij from the chargino-sneutrino loop is given by 

^'^'''' = -^fA:V), (7) 

whereas the Z-contribution reads 

^ IGir'^g'^ sin^ Ow^'z ^ ^ ^ ^ ^ 

with X = {LL, LR, RL, RR}. In the above, Oy'^ denotes combinations of rotation matrices 
and coupling constants and s(x'^) and t(x^) represent the Passarino-Veltman loop functions 
which depend on = m|_/m? (see [IS]). Since <C m?, the Z-mediated contribution 
would dominate over the one coming from 7-penguins. However, in the MSSM, this is prevented 
by a subtle cancellation between different diagrams contributing to the leading Z-mediated 
contribution [15]. It can be understood by focusing on the dominant contribution to ii 3ij 
which is given by diagrams where the external leptons are left-handed (the other cases are 
suppressed by the Yukawa couplings of the charged leptons). It is given in [lOj by the form 
factor 

Fll = , (9) 

g"^ sin^ 9yi/m?2: ' 

where zf' = —-^{—\ + s^^) is the Z — ii — ij tree-level coupling (i = j) and is the Z — li — tj 
1-loop effective vertex, with i ^ j. In the MSSM, the sneutrino mass matrix is diagonalized by 
Zy, a 3 X 3 unitary matrix. But when chargino mixing is neglected, the masses cancel out in 
the combination of loop functions from different diagrams and one ends up with F'^ oc {ZyZv)ij 
which vanishes for i ^ j due to the unitarity of the Zy matrix [15j. However, this cancellation 
can be spoiled if new interactions are introduced like it is the case with the inverse seesaw 
mechanism-|-. In particular, the — i^R — coupling gives rise to new diagrams containing a 
higgsino-sneutrino loop which spoils the cancellation since they give Ff^ oc Zy Zy*Y^^*Yi!^-' . 
It is worth noticing that this does not happen in the MSSM extended by a high-scale type-I 
seesaw since the very heavy right-handed sneutrinos suppress the higgsino-sneutrino loop. In 
the super symmetric inverse seesaw, the large neutrino Yukawa couplings and the TeV-scale 
right-handed sneutrinos make this contribution naturally large. 

The complete derivation of the analytical expressions for can be found in [lOj and only 
an abridged version will be presented here. The chargino-sneutrino loops can be decomposed as 

i7(0) L_ f + ^ f 10) 

Considering that left-right mixing is negligible in the sneutrino sector when the neutrino 
Yukawa couplings are relatively small (which will be verified by our numerical study), the two 
contributions reduce to 
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Table 1. Current experimental bounds and future sensitivities for different cLFV observables. 



^i°wl„o = jh,{3gcw+g'sw) = ^Sii(l + 2cl,), (11) 

It can be noted that the wino contribution vanishes for i ^ j, which corresponds to the 
cancehation discussed earher, while the higgsino contribution, which is absent from the MSSM, 
does not vanish, explaining the different behaviours between the two models. We have verified 
that the expressions above reproduce the numerical results based on complete, non-approximate 
expressions. We have also verified that the higher-order terms in chargino mixing are numerically 
negligible. It is worth mentioning that the F^^ vertex does not depend on supersymmetric 
masses. As a consequence, the Z-penguin exhibits a non-decoupling behaviour and a large 
SUSY scale does not suppress its contribution to cLFV processes. 

Finally, although the previous discussion has been focused on ii — )• 3ij processes, the same 
enhancement in the Z — ii — ij effective vertex also affects other observables. This is the case for 
fi — e conversion in nuclei [T6] and the case of r ^ P^ii, where P° is a pseudoscalar meson [17J. 



3. Numerical results 

In this study, we focused on leptonic observables which can be divided in three broad classes: 
radiative decays, like /i 67, 3-body decays, like r — > 3//, and neutrinoless conversion in muonic 
atoms, like fi, Ti — )■ e, Ti. These observables are actively searched for at many current and future 
experiments. The list of processes we considered together with the current upper limit and 
expected future sensitivity can be found in table [TJ We also addressed high-energy observables 
such as flavour-violating neutralino decays Br(x2 ~^ Xi^i^j) ^^'^ slepton mass splittings Am^ 
(m^- — m^- ). However, we found that the small neutrino Yukawa couplings needed to comply 
with experimental limits on /i, N 6, N conversion (as discussed below) gives low rates for 
X2 ~^ Xi^i^j s-iid small mass splittings. We have also verified that cLFV decays of the Z and /i° 
bosons are not significantly enhanced. 

Our numerical results have been obtained with a SPheno |29[l30j code generated with the 
Mathematica package SARAH |31H33j . The computation of the cLFV observables is based on the 
results presented in [T3l[16l|T7] , that we extended to the inverse seesaw case. As such, it is not 
limited to the Z-penguins but also includes contributions from Higgs-penguins, 7-penguins and 
supersymmetric boxes. The analytical expressions for fi — e conversion in nuclei, Br(£j — )• 3ij) 
and Br(r P^ii) in the case of Z-penguin dominance are given in [1(}]. The supersymmetric 




Figure 1. /x — e conversion rates in ll'Ti (left) and Br(/x — > 67) (right), as a function of M 
and iYjY^)i2 for different M/? values: Mr = 100 GeV (blue), Mr = 1 TeV (red) and Mr = 10 
TeV (black). We set Aq = —300 GeV, tan^ = 10, sign(;u) = + and Bq = 0, and we randomly 
vary mo and M1/2 ™ the range [0,3 TeV]. The horizontal dashed lines represent the current 
experimental bounds and the dotted ones represent the expected future sensitivities. 



parameters of the model are niQ, M1/2, AQ,tan f3, sign(//) and Bq. Since they are of little relevance 
for Z-mediated observables, we fixed sign(/i) = + and i?o = in the numerical analysis. As 



can be seen in figure 1(a) , when an observable receives its dominant contribution from Z- 



penguins, it has a very small dependence on the cMSSM parameters making our results quite 
general. However, this is not true when another contribution is dominant like in Br(/x — > 67), 



see figure 1(b) for example, in accordance with the literature [21 The results presented 
here are obtained with degenerate singlets Mr = diag{MR, Mr, Mr) and we have verified that 
processes dominated by Z-penguins are only marginally affected by non-degenerate singlets. For 
example, the relative change in the conversion rate for ^u, Au — )■ e, Au is always below 2%. As a 
consequence, the results discussed here hold true for both degenerate and non-degenerate right- 
handed (s)neutrinos. Since fix does not explicitely appear in analytical expressions for cLFV 

processes, the relevant parameters are M ~ (as defined in Eq. ([5])), which controls the size 
of the Yukawa couplings, and Mr which controls the size of the non-SUSY contributions. 

In figure 2(a), we display the absolute contributions to Br(/i — t- 3e) in order to illustrate 
the non-decoupling behaviour of the Z-penguin as previously discussed. We define the relative 
contribution to cLFV observables as 



(13) 



where Rj is the rate (branching ratio in case of ^ 3e and conversion rate in case of // — e 
conversion in Au) that would be obtained, should the i-contribution to the process be the only 
one. It can be seen from figure 2(b) that Br(^ — > 3e) is clearly dominated by the Z-boson 



contribution in most of the parameter space, the only exception being for very low niQ = M1/2 
where the 7-penguins dominate. This can easily be understood by recalling the scaling of the 
operators for the different contributions in equations ([7]) and ([8]). The same behaviour was 
observed in other observables as described in [10] and allows to conclude that, if present, the 
Z-mediated contribution will dominate cLFV observables in the super symmetric inverse seesaw, 
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(a) ^ — > 3e (b) fj, — e conversion in T-g^Au 

Figure 2. Absolute contributions to Br(/x — > 3e) (left-hand side) and relative contributions to 
fj, — e conversion in 79^ Au (right-hand side) as a function of mo = foi' & degenerate singlet 
spectrum with Mji = 10 TeV and M = 10^^ GeV. The rest of the cMSSM parameters are set 
to ^0 = —300 GeV, Bq = 0, tan/3 = 10 and sign(/x) = +. On the left-hand side, solid lines 
represent individual contributions, 7 (black), Z (blue) and h (red) whereas the dashed lines 
represent interference terms, 7 — Z (green), 7-/1 (purple) and Z — h (orange). Note that in this 
case h includes both Higgs and box contributions. On the right-hand side interference terms are 
not shown to make the results clearer. The individual contributions are 7 (black), Z (blue) and 
h (red) and boxes (purple). 



except for a very low SUSY scale. 

When the Z-penguin gives the dominant contribution, there is a strong correlation between M 



and the processes rate. This can clearly be seen in figure 1(a) for the fj, — e conversion rate. Prom 



the non-observation of cLFV, an upper bound on M, or equivalently on the size of {YjY,y)ij, can 
be derived. Among the many observables involving fi — e transitions, the most constraining one 
is neutrinoless conversion in muonic atoms and the approximate limits it gives are summarised 
in table [2] for a conservative choice of Mr > 1 TeV. Smaller values for Mr will only result in 
stronger limits, especially when non-SUSY contributions become dominant, making — > 67 the 
most constraining observable. It should also be noted that future experiments like COMET, 
PRISM/PRIME and Mu2e have the potential to greatly improve these limits. 
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Table 2. Approximate limits on M and {YjY,y)i2 from the non-observation of /x — e 
conversion in nuclei. This table includes current experimental bounds and future expected 
sensitivities 



Other cLFV transitions also exhibit a clear enhancement due to the Z-penguin dominance. 
In particular, the branching ratios of different ii — )• 3ij channels are related by 
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W^2 _5,^^. 

(14) 
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From this equation, we can expect that, considering similar combination {YyYy)ij, all the 
branching ratios for £j — > processes will lie with 1 — 2 orders of magnitude. However, 
the experimental upper limits on these observables are more relaxed for r — e and t — /U cLFV 
transitions than for fj, — e cLFV transition. This implies that 3-body r cLFV decays have little 
possibility to be observed in the near future, except if there is a strong cancellation in Y^Y^, 
which would suppress fi — e cLFV processes. The same conclusion also holds for semileptonic r 
decays. Finally, it is worth mentioning that in Z-penguin dominated scenarios, there is a clear 
correlation between the rates for ^ 3e and jjL — e conversion in nuclei. Numerically, we found 
CR(// - e, Ti) /Br ^ 3e) ~ 15. 

Z-dominated observables are quite insensitive to the cMSSM parameters, implying that 
these results should hold when the recent mass constraints on the Higgs boson are imposed 
on the model. We have checked this using a benchmark point associated to a Higgs mass 
of m/jO = 126.5 GeV and the corresponding predictions for the different observables are given 
in [To]. Moreover, since the neutrino Yukawa couplings are constrained to be below 10~^, we 
do not expect the singlet sector to sensibly modify the Higgs mass and renormalisation group 
equations (RGEs). 

To conclude, the super symmetric inverse seesaw is a very attractive extension of the MSSM 
since it generates neutrino masses while keeping all the new Physics below the TeV scale with 
naturally large Yukawa couplings. In this work, we considered its predictions for different 
lepton flavour violating observables mediated by the exchange of a Z-boson and found that 
the Z-mediated contribution dominates in most of the parameter space, providing an important 
enhancement to different observables. As a consequence, the most constraining observable is 
now neutrinoless /i — e conversion in nuclei, in contrast with the usual behaviour of the MSSM 
extended by a high-scale type-I seesaw. This allows us to put constraints on the size of the 
Yukawa couplings {YuYy)i2 > 2.7 x 10~^, clearly excluding Y^ ~ 0{1) as usually considered. 
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